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Abstract For the wave representing particle traveling through any layer system we calculate 
appropriate phase shifts comparing two methods. One bases on the standard scattering theory 
and is well known another uses unimodular but not unitary M-monodromy matrix. Both meth- 
ods are not equivalent due to different boundary condition - in the one barrier case there exist 
analytical expressions showing difference. Authors generalize results to many barrier (layer) 
system. Instead of speaking about superluminarity we introduce into the quantum mechanics 
so called by us "hurdling problem": can a quantum hurdler in one dimension be faster then a 
sprinter (without obstacles) at the same distance. Relations between wavefunction arguments 
and delay or advance are shown for Nimtz systems. 

1 The tunneling times definitions with reference to S (scatter- 
ing) and M (monodromy transfer) matrix theories 

1.1 Smith's method as S-matrix method 

Before 1960 duration of a collision was a rather ill-defined concept, depending on a more or less 
arbitrary choice of a collision distance r. Such a point of view was represented by F.T.Smith 
(1960)[1] in his paper "Lifetime Matrix in Collision Theory". In that work the author tried to 
generalize delay-time At = h(drj/dE) resulting from analysing the scattering of the wave packet 
into a concept of the general S matrix theory according to papers written by Bohm (1951) [2] and 
Wigner (1955) [3]. If collision time is defined as a limit for r — > oo, then the difference between 
the time: 

a) in which the interacting particle stays within distance r, 
and the time: 

b) it would have spent there in the absence of the interaction emerges as a well-defined quantity 
which is finite if the interaction vanishes rapidly enough at large distances r ^ oo . 

"In quantum mechanics, using steady-state wave functions, average time of residence in the 
scattering region is the integrated (excess) density divided by the total in-or out (ward) flux, and 
lifetime (more precisely, time delay) is defined as the difference between these residence times 
with and without interaction." 




average integrated density / flux 




where average density in the absence of the potential is: 

Poo =< P(.x) >=< i^*ooix)'4'oo{x) >= lim - / 
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and jinw/outw is the inward or outward flux as defined by Smith 

fik 

jinw — = joutw = (1-3) 

m 

where, according to the scattering theory, the asymptotic (one-dimensional) form of ip at large 
X is: 

i;o^{x) = A{e-''''' - e'^'^e'''^) (1.4) 

and A is normalization while density in the central region is: p{x) = ip*{x)'ip{x) 

ipoo in case of many channels and separable radial part of the Schrodinger equation can be 

written as 
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FIGURE 1. Inward, outward, tunneling etc. waves in scattering. During scattering only one wave 
(here) outward or inward is modified. There is no cross terms between inward and outward fluxes . 
'^tun is not incorporated in vj/o"*"" nor in It is not clear if the reflected wave ^refii- \S\ = 1) 

is equal to \I'°"*'". In shadow region the complete wave function must vanish [16] {ip = 0), there 
is no place for "if tun in V'- 

If the wave functions are normalized to inward or outward unit flux through a sphere with radius 
r — > oo,than on the basis of complete wave functions (cf.fig 1) we build the lifetime matrix Q, 
using the time operator (there are no consistent theory till now concerning the time operator ) 
t = -ihd/dE 



Q 



(1.6) 



where S is the scattering matrix. According to Smith's paper, Q and S contain complementary 
information and after diagonalization of Q its eigenvalues are the lifetimes of metastable states, 
while the corresponding eigenfunctions are the proper functions describing these metastable 
states. That's why Q is called the lifetime matrix according to the formula derived by Smith as 
below: 
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where the average value is taken to ehminate oscillating terms at large r. Q is introduced 

corollary using identity Q = Hdrj/dE. 

Ohmura generalized above consideration on time packets : 



sc 
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In his method A,a,/3 are real functions, da/dcu gives time delay of incoming time packet due to 
reshaping before and during collision while d{[3)/doj due to reshaping only during collision (/^ 
is the differential cross section). Using time dependent flux formula j{t) averaged over time: 



The above idea has been applied by Olkhovsky-Racami [4] in investigations of reflection and 
tunneling times. All these methods analyze variations of the complex wave arguments during 
scattering, directly (by Ohmura) indirectly in terms of fluxes in ref.[4]. Below we try to flnd 
analogue of phase functions df3/du (distributions) as function of projectile wave-number for 

transmission through systems as in Nimtz superluminar experiments [12]. 

Depending on the problem under consideration the scattered phase shifts can be defined in 
reference to other known shifts (as Coulomb phase shifts or just kr the argument of undistorted 
i^mw(^outw -^aves i.e. - without interaction phase at r is given simply hy kr). Now having the 
scattcrcr we replace it by potential (repulsive or attractive) and matching wave functions and 
their derivatives outside potential range (as solutions of corresponding wave equation with initial 
condition that the wave function is equal to zero at origin as well with assumption that both 
fluxes inward and outward are orthogonal [1], cf. figl) we calculate scattering amplitudes. The 
scattering device together with incoming flux is located at cm. and elastic channel is one usually 
created by the nonresonant "reflected" wave function with the same k vector. The scattering 
theory doesn't make difference between elastic reflected and transmitted waves. There is only 
one averaged elastic channel wave function. Such situation is typical for all 5-matrix problems 
in area of nuclear reaction, the phase shifts define scattering amplitudes and these quantities 
define cross sections to be considered. The phase shifts are not monotonic functions of energy 
[5] and such dependence were not investigate due to not unique definition of potential. There 
were attempts to solve the inverse scattering problem (from phase shifts to restore potential) 
but without success. 

The incoming flux when scattered by the target (barrier) is converted into the outgoing parts i.e. 
reflected and transmitted. On the projectile side in one dimension thought experiment there is 
reflected particle interfering with incident beam while on the other transmitted. But in reality 
it is not easy to say which particle is reflected or not. In the stationary theory we take into 
account only an averaged outgoing flux (mixture of reflected and transmitted particle; cf. figl, 2.). 
From S- matrix point of view we have in one dimension two subchannels (R,T) or as in case 
of the separable radial part of the Schrodinger equation we must remove 1-wave degeneration. 
The 1- wave splits into two subfunctions corresponding to the reflected and transmitted 1-wave 
functions. In case of the radial coordinate (one dimension in three dimensional space ) we are 
unable to define the left - right sides even experimentally to distinguish the reflected wave from 
the transmitted one (the exception is the shadow region in flg. 1). This degeneracy cannot be 
describe with traditional Si elements. The scattering on the set of many barriers treated as one 
"black box" should be described by S'-matrix. But such system should be characterized by one 
phase shift rj or rji what cannot be true. We have two functions in output each with its own 
phase shift. To describe such system we must introduce unimodular M- matrix. The M-matrix 




(1.9) 



he got the mean time delay: 



Ai = 



jA^P{u)vdu JA^vdio 



(1.10) 



conserves the mutual exclusion relation between amplitudes R and T. From the transfer matrix 
point of view each wave function should receive proper phase shift {ipR^ipT ) after scattering. 
Then we can find S- matrix amplitudes if M- called transfer matrix is known. There is one more 
problem Si phase shifts are found from one matching while in case of transfer matrix, it is built 
from multimatching conditions as consequence of many borders between media (inside nuclear 
structure). The scattering system is no more the "black box" type. The internal structure causes 
multiscattering as sequence of subsequent reflections and refractions what requires description 
in terms of not unitary M- matrix. The question is if both descriptions in terms of ( S" if exist 
and M) are equivalent. There is no papers where S- matrix R,T amplitudes could be calculated 
independent of M. (5 is deducted from M not vice versa). 

Till now the transmission (tunneling) was taken into account indirectly through reaction chan- 
nels. Such treatment put tunneling outside the scattering theory. The consequences are : 

- Non unique solutions at the origin (from two solutions we take only regular one into consid- 
eration, analysis below shows that in tunneling case the wave function is different from zero at 
system origin or at least undefined-not used) ^. 

- In nuclear physics there are problems with hard or soft core potentials which were not tested 
or compared with tunneling. 

- Validity of time reversal invariance or detailed balance theorem which says that the time 
reversed incoming state (under the operator K) is equal to an outgoing state with the same 
energy. The reversed in-state goes into the asymptotic free time reversed state when t — > +00 
i.e. k'^f'^ = ^'ji"-' and = *jt^ These relations induce S'j/// = 5/j called detailed 
balancing or microreversibility. In other words the transition probability for the inverse process 
with time-reversed parameters is the same as that of the direct process. 

But tunneling is irreversible process and we suspect- cannot be described by function regular at 
origin. 

In general case of the reaction a + A h + B (in the subbarrier collision) tunneling in out-state 
{h+B) is different from that in in-state [a + A) and tunneling disturbs scattering states. Clearly 
nonunitary condition breaking microreversibility relation , tunneling however can be introduced 
as additional indeterminance in scattering theory. 



1.2 Weak Wigner causality and Wigner time 

The Wigner time is the simplest one. According to [3] and formulated there the principle of 
causality, the scattered wave cannot leave the scatterer (of diameter r) before the incident wave 
has reach it i.e. dr)/dk > —r. This expression in case of positive derivatives gives retardation 
while negative values advanced in time solution, for the outgoing wave as defined in [3] we can 
write tout = ;jj + §§f • Experimentally it is not easy to find from the excitation functions (cross 
sections) r7(fc).(cf.eq. 1.10). In reahty in macro world the scatterer (Coulomb or gravity field) 
has infinite radius what forces tout 00 • Let d = 2r and tin = then 

_ _ _d 2dr] 

twigner — tout tin — ~^ vdk \ ) 

If we know r]{k) the Wigner time (the group one) can be easy derived for finite systems. 



1 . 3 Mono dromy 

To introduce M matrix we need two ingoing m(ib) and two outgoing out{±) particle wave 
functions, (cf.Fig2). 

FIGURE 2. The monodromy problem as defined in [6,10] for one barrier. Mutual relations 
between ingoing and outgoing (from left or right side) particle wave functions are displayed. 
(In the picture bars mean complex conjugations.) The transition from initial state to final 

^Jost functions start from two irregular solutions but physical meaning has their sum equal to zero at origin 
[5, ch.ll, eq.(5,26,71,72)] i.e. sum is regular . 



one ( I ^ ( I is given by the unimodular M matrix I | easy 

deducted from transitions as drawn in picture. 
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The transmission or reflection through any periodic or aperiodical set of square barriers rewrit- 
ten as the transformation from to (<&t) (undergoing the monodromy matrix), can be 
described in terms of the [2 x 2] transfer ("monodromy") cells given by superposition of [O^] 
matrices with [Hi] . [Of] represent free wave propagation between barriers which can be in- 
terpreted as phase translation to given position (the middle and/or edge of barriers) or phase 
translation about relative distance if [O^^j^] are used . In ref. [7] authors use simply the name 
translation operator. [Hi] describe particle motion under any barrier and are responsible for 
interactive wave propagation. Superposition of both matrices [Hi] [Oi] represents propagation in 
two opposite directions inside or outside media. In case of the square barriers unimodular [Hi] 

matrices have form ^ where a, b, c are real, 
c a 

Monodromy (or the transfer matrix or translation operator unimodular not unitary (as well 
not equivalent to unitary) transforms the initial wave function amplitudes Ao,Bq or (1, i?) onto 
outgoing one An,Bn or (T, 0). 
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The monodromy form of M depends on the basic wave functions to be chosen. M shifts the 

solution of the Schrodinger equation from x to x + d i.e. from beginning of the barrier system 
to its end. In the time depended approach the wave function underlies unitary evolution : 
^(t = -|-oo) = U{+oo, —oo)^in{—oo). The initial and final wave functions are separated in time 
in S-matrix treatment while in M matrix approach as well in space : far left, far right. 
For inward, outward [(complex exponential functions, Hankel functions etc) or real basic solution 
like (cos, sin, regular and irregular Coulomb or Bessel functions etc.) , here both representation 
M and M' are connected by the unitary transformation] we get: 



(1.12) 



Monodromy M represents propagation of the wave functions through system of multiple cells 
as sequence of reflections and transmission (at each cell edge the wave is splitted into reflected 
and refracted (tunneled)) or M can be interpreted as superposition of cells characterized by 
two waves inward and outward. The four M- matrix elements can be expressed as function of 
complex variables T and R, above relations define only M21 = R/T , M22 = 1/7" elements. The 
remaining M12 , Mn elements, connected by det[M] = 1 relation, we deduce from matching 
conditions. The monodromy is unimodular not unitary. 7^ and hermitian conjugation 
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does not describe inverse motion. Multiple reflections and transmissions are strictly correlated 
with multiple matching. If system is asymmetric (there exist at least one left and right matching 
which do not coincide), equivalent S'-matrix can exist if we introduce additional phase shift A(/p 
between R and T waves. Monodromy can be periodic but not necessary. Without dissipation 
(energy loss) system consist of multiple superposition of unimodular matrices. 
Now we consider transmission through certain device created by superposition of many barriers. 
Such systems can be equivalent any arbitrary shape potential U{x) defined on the intervals 
— €-1 < X < Qi + ei with help of square barriers (e.g. barriers on the Cantor set etc.). There 
are barriers as in the fig.3: 



-2£- 



a. a. a^ distance 

(barrier positions) 

FIGURE 3. Set of barriers as used in tunneling and transmission or reflection. 
The superposition of all "phase translations" describe the full transfer operator M as transfor- 
mation from the initial (spinor) amplitude state to final one by means of matching conditions. 
By appropriate unitary transformation we can choose convenient amplitude representation. The 

choice depends on physics to be considered. 

The n-th barriers system can be described by barrier center coordinates (oj) or interbarrier 
distances (Aj) and barrier widths and heights {2ej,K^p). Then the system width is given as 



(1.13) 



ci = e„ + cin - 0-1 + £i or d = J2j {2£j + Aj) 
O'k+i = + fflfe = J2j dj 

where is the i-barrier center position and 2ei width of the i-th barrier, Aa^ - interbarrier 
distance (between barrier centers) and Aj free cell width (distance between neighbour barrier 
edges. 
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matrices. The 



The M transfer matrix can be expanded as multiplication of [-ffjjand 
first discribe particle motion under the barrier (or in media) while second free motion between 
barriers. The matrices [0]j can be written as the function of Aj i.e. distances between adjacent 
barriers given by the difference of their edge positions: 



(1.14) 



The transfer matrix can be expressed in terms of the barriers edge i.e. it Cj or distances 
between adjacent barriers eq.(1.13) then the transfer operator M' is 



n-l 



m = K] n [oaj [Hi, 



(1.15) 



i=l 



In representation such that [OaJ is diagonal, we have the final form of the position independent 
complex transfer operator 

Um [M'] Ul = [M] 



in case of the cos, sin base U is the unitary matrix which make diagonal [OaJ i-e 
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is interpreted as stream of two waves propagating in opposite directions. 



All that matrices, denoted by 



to monodromy type [6] ,(cf fig2.) i.e. 



, like 
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, [M], and [O] after diagonalization, belong 



X + iY -V-iW 
-V + iW X-iY 



(1.16) 



If [M'] is real and [M] = [M'] C/]^ then X = {M[^ + M^2)/2; W = {M[^ - M^2)/2; V = 
(M(2 + M^i)/2; Y = (M{2 - M^i)/2; 

V and W gives information about asymmetry in tunneling (breakdown of balance between two 
waves traveling inside the "black box" in two opposite directions). 

The product of O^'' H^^^ = Mi is then the element of the barrier structure named the single 



cell transfer operator. 
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t can be written as: 
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^ + CTj) sinh(2Kiej 



cosh(2Kjei) + f (^ - (Ji) s'mh{2Ki£i J ^Va, 

+ "^i) sinh(2Kiej) cosh(2Kiej) - |(^ - CTi) sinh(2Kiej 



(1.17) 



where = and = {nf^)'^ — fc^, k is the projectile momentum. Transition to 5 barrier set 
is done using eq.(1.17) when in Mj we put A = 2K|e then we get: 
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using two phase representation (of reflection and tunneling 
for symmetric structure in more compact form: 



we can rewrite the [Mj] matrix 
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(1.18) 



where we put on the base of the single barrier transmission formula : 



'^Monodromy, i = + 9^2,1 = tan {^{— - (Ji) tanh(2Kiej)) 



(1.19) 



reseparation of both phases in general case (many barrier system) is not easy. 

This single cell operator can also be written with help of amplitudes (-Ri,i,ri^j) and phases 

(</?2,i,¥'i,j) as: 



Mi = 



O 



m 



JkAi 







-ikAi 



1 p(«'/'2,i+¥'l,i) 



Tl 



(1.18) 



We assume Ti = Ti,ie^('^2,<+¥'i.i) ^nd Ri = Ri^ie'^f2,i+'Pi,i+^'Pi) 2 go the wave phase in the reflec- 
tion channel differs by A(fi from transmitted one. Acpi can be computed from the expression 

"^^^ (real 



Due to relation |_R| + [r| =1 all cyclic function can be calculated modulo tt i.e Ti,; 



ITi Je^ 



(cos) sinusoidal amplitude= ± modulus) while the phase between R and T amplitudes due to imaginary factor 
i modulo 7r/2 (equivalence of tan(na:) and cot(na;) sets, see [2] cli.ll, Bolinis ipi are different from ours, his 
Vrefl = ^trns ± 7r/2) ; We require only smooth behavior of phase function and its derivative. 



tan(A(/?) = y 

Aif ^ for asymmetric systems, for symmetric one we can put A(p = 



(1.20) 



1.4 Phase properties of M matrix-one barrier case or the barrier set equiv- 
alent to one 

Let in the equation (1.17,18) A = constfy = e and H is the same for all i Hi = iJg. Then we 
call Mi = Mg^A single cell "power" (or periodic) monodromy operator. Internal structure of [M] 
representing certain device causing reflection and refraction defines transmission or tunnelling 
through the barriers as well as general (aperiodic) monodromy. We maintain that monodromy 
as applied to two channel elastic scattering is group property . 

Let we make one more comment: periodical structure emerge with multiple application of 
[Oa] but final boundary condition can change periodicity. 
We can solve eq.(1.12) with M as in (1.16,17) to find amplitudes R and T 



= Msiieo 
The general solution of that equation is: 
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(1.21) 
(1.22) 



We can calculate also the amplitude ratio 



R 
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In above formulas we put 

M22 = X - lY 
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-f exp(z(A(^ + k{an + + ^i))) 
-'1 



^ exp(-z(yi + (^2)) ^ (iX.p{-iip Monodromy) 



and system total width d = a„ + — ai + ei can be expressed by intcrbarricr spacing Aj 
i.e.: d = J^iC^^i + ^i) = J^idi- In that way di defines the single cell width (one barrier plus 
one interbarrier well). As we have seen in explicit formulas for and M22 the argument of 
diagonal matrix elements is arg(Mj^) = <fi + <f2 = ^Monodromy ^"^^ it depends through ipi on 
internal system structure. 

Usually for symmetric systems we have 1^ = then V is real and both amplitudes have the 
same phase. 

The expression (1.21,22) for amplitudes R,T depend on the wave function value at the "black 

box" edges. In case of symmetrical aperiodical systems we put fln + = + ^i- 
For one barrier using unimodularity, in case of R and T we get: 



iT 



e ^'^'^ = zsin((^i) exp z(tan ^ — ai) tanh(2«:i£i) — fed) = isin((^i) exp(i(^2) 
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(1.21) 
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cos((/?i) exp ?(tan ^ — ci) tanh(2«;i£i) — kd) = cos((^i) exp{i(p2) 

(1.22) 



We can rewrite i~p2 as follows 



= + 2tan ^{ — tanh(Kiei)) + tan ^— . = 2ry + (1-23) 

o\ (1 + erf ) sinn(2K;iei) 

Here r/ is as in [9] (Aufgabe 57) i.e. r\ = —kd/2 + tan^^ tanh(Kiei)) , Fliigge to solve 
problem put V'(O) = 0) '^^ do not need that condition and "translational" boundary conditions 
at r = ro =b d/2 (assuming rg = 0) result in additional phase ^pi. Then in our methods cross 
section is proportional to sin^((^2/2) not sin^(r;) . 

In fig. 4,5 we have shown transmission , 5 argT = (p2{k) and 9^2 (^) = 2rj + kd in single barrier 
case. The ^2{k) in allowed k-band is increasing function of k and the quantum hurdler is slower 
then particle without obstacle; the r]{k) function has not such properties. 

Using ipMonodromy — Et = 0, from position of the packet center we find the transmission time in 
the Nimtz experiments. 

The heuristic time calculation fulfills typical limits : 

fid{ifi^i + ip2,i) ^ 2_ 

idE KE-^oo yn- 

For typical ip2 expressions we get the same limit. That expressions are not here important. We 
presume wide systems are composed from thin elementary segments. 

FIGURE 4. Transmission through one barrier three units [mm] wide, barrier height in k 

units [1/mm]. —kd represents maximal negative phase slope according to weak Wigner causality. 
Slope of (p2 = V'2 + ~ is connected with the group velocity in transmission through the 
barrier, {d / vdk){(p2) = 5t is time "delay". The wave length A at k^'^-* is smaller than the barrier 
width. 
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FIGURE 5. Transmission through one Nimtz barrier 6mm wide, k*^*^-* is barrier height in k 
units (here mm~^). As in fig. 4 (—kd) represents maximal negative phase slope according to 
weak Wigner causality. The wave length A at is bigger than the barrier width so phase 
characteristic is dominated by ipi. The slopes of (pi{k) and 5 argT phase curves are positive and 
give retardation. However ip2 — kd = 2r] suggests small speed advance. 
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In Fig. 4 we have shown the phase characteristics in tunneling through one (or two barriers see 
fig. 6). Due to weak Wigner causahty apphed to the sum of both phases {d / dk){ip2) > —d. From 
monodromy 6 arg T = 5 arg R = + ifi — kd = (p2- The phase ip2 — kd alone is typical 2r/o as in 
Q.M. -textbooks (see [5,9]). Sometimes for -one or few barriers - when A is comparable with the 
barrier width, the height of the barriers can be easy deduced from the phase characteristics. It 
is not a rule cf. fig 5, 7 and others. 
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FIGURE 6. Transmission through two equal barriers each one length unit [mm] wide. The 
cavity diameter is one millimeter wide too. is barrier height in k units [1/mm] the same as 
for one barrier tunneling (cf. fig. 4). 

We can write general expression for monodromy single cell traces 



cos (pi = \Tr [Mi] = Re 
= cot(v9i ,) cos{^2a) - sm((/?2,i) 



exp(-j((^-^_.+ip2 .)) 



cos(y-^ .) 



(1.24) 



For each cell we can define two internal phases (^i^j, (p2,i "Bloch phases" ( in analogy to (pi 
in [10]) and another one /cAj as external typical for the interbarrier movement. Such method 
can be compared to the scattering as in [10]. Stability in classical mechanic is expressed by 
inequality TrM < 2, here in quantum mechanics TrM describes mutual ratio of refiection and 
transmission (tunnelling) - expressed respectively by (pi or ^pi behavior. High above the barriers 
^1,1 ^ /"^ so (pi ^ ip2,i — > kdi- We see that much more appropriate phases to be named "Bloch 
phases" are (pi^i,(p2,i which describe internal device structure however here we deal with the 
aperiodical or quasiperiodical system (periodic H-boundary). 

Let assume there exist average transfer operator M equivalent to superposition of equal or 
different elementary monodromy cells (assuming symmetric case an + Cn = — oi + ei)- We say 
there exist an equivalent "black box" barrier operator M which preserves the single cell form 



(cf. eq. (1.15,17,18)). 
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(1-25) 

The phases (pi,(p2 are functions of device internal structure. We assume such phases exist and 
can be computed while M is folded from square barriers, Te*^^ = f^^{d,Q = sm(ipi)e'^'^^— 
where C represents all internal variables needed to compute ^1,^2 ■ most cases the product 
of the single cell expressions (1.18) can be computed only numerically. The total phase shift 
change is given approximately by (M transforms input (1,-R) onto (T, 0) output): 

^Monodromy = + ^) = arg{/(d C)} = " tan~\^{^ - a) 1^11^2^.))'' 
Equivalent "black box" width is d = 2e . Well above the barriers i.e. arg{/((i, Ci'^))} ^ 

a^i=\/—l 

there is the transmitted wave running in the initial direction nearly without distortion. Here 
the argument of tan~^ has only symbolic meaning however we assume existence of a,e,K for the 
equivalent black box barrier. 

In equations (1.25), for periodic case (i.e. Mj = M^^a ) where Cayley Hamilton theorem can be 



applied. 



O 



([/)■ 



just cancel exponential term on the diagonal of Mj matrix. The total transfer 



periodic matrix is if cos((?!)) = ^TrM^^A- 



M 



O 



-1 



A 

sin(Af(/)) 
sin<i 



<A 



o 



M. 



sA' 



sin(JV(/>) 
sin<^ 



^ sin((jV-l)0) 
sin</> 



o 



sin((Af-l)(^) 
sind) 



it is obvious that similarity of ^ trace with M trace (as in [10]) is accidentally. The trace 
properties are useful when we want to write any power of in terms of Mg a and the unit 



matrix (see ref [13]). There exists however an additional matrix factor 



O 



which changes 



the final trace completely. From physical point of view we are interested only in transformations 



which put 



O 



into diagonal form meaning that between interaction areas we have two free 
waves running in opposite directions. In general case for the i-cel\ we can define two scattering 
or Bloch phases ip2,i, <^i,i (Bloch phases suggest periodicity what is not here the case) and phase 
displacement kAi so each cell has different trace properties. During out of resonance tunneling, 

particle seems to be insensible to Aj distances. 



It is impossible to make diagonal both 



O 



and 



(U) 



Matching conditions induce the 



monodromy form of 



(U) 



describing particle movement under the barrier as in [6] . The total 
transfer matrix is composed of many matrices. We must know its final form explicitly to find 
mutual ratio of reflection and transmission. 

In aperiodical symmetric system (like barriers on the Cantor set) the averaged term ipi + ip2 
is different from the expression valid for the single barrier (1.21,22) . The case symmetric is 
important, it is easy to calculate the phases ipi,(p2 '■ in the monodromy matrix (1.16) W = 
and as written already V is real so from (1.24) V = cot((^^) = R/T. Prom X, Y we extract the 
second phase. For asymmetry we get cot{(p^) = ±\/^^ + = R/T and tan(A(^) = A(p{k) 
is phase difference between amplitudes in reflection and transmission. 
Resume: 

Translation operators suggest "translation in time" too. S-matrix is time independent opera- 
tor {U{t — oo,+oo)) and treats the quantum wire translation device as "black box" which 
structure should be find out in "phase shifts" experiments. Unitarity of S'-matrix suggests full 
symmetry of "black box". It is not clear if S'-matrix phase shifts r]i can be used to calculate 
time delay for particle traveling through the investigated object cf.[l]. 



Only the transmitted (tunneling) waves "feel" the object size (i.e. its depth). Tunneling intro- 
duces asymmetry into experiment as well into theory. May be motion emerges in the quantum 
mechanics as consequence of reflection - transmission interference (which takes place only under 
/ over the barrier or the potential well). So there exist equivalent iS-matrix {Sm) related with 
M-monodromy translation operator if M is symmetric. Asymmetry causes problems and mi- 
croreversibility is exact up to tunneling. M-matrix can be periodic but generally it is not. To 
describe properly tunneling and reflection we need mixture of inward and outward solutions (or 
at each point of the space the true solution is mixture of regular and irregular one coupled by 
tunneling effect at origin where integration of the wave equation starts). 

In tunneling, due to equivalence between the complex Schrodinger and two dimensional Maxwell 
(Helmholtz) equations, we consider as well photons as massive particles. 



1.5 Monodromy time 

Two phase shifts suggest that Wigner causality should be applied separately to reflected and 
transmitted waves. But if reflection has nothing common with transmission their coupling 
through the |i?p-|-|rp = 1 relation would result in completely different phases of both functions. 
The above analysis shows that in principle for symmetric systems both waves have common 
phase (p2 = ^1 + ^2 — kd [8] and such sum should be used in causality relation. The phase ip2 

alone corresponds to rj (rji) if k ^ k^^^ where k^^^ is the barrier height i.e. ^2mV/h^. We can 
generalize the Wigner time and write 

tMonodromy = ^^^^dk''^ > 0; <P2 = + 27/ (1.26) 

Tangent to f Monodromy = 'f2{k) + dk cannot be negative: the scattered wave cannot leave the 
barrier of width d before the incoming wave has entered it - in consequence tMonodromy > 0. 
There is additional phase (pi which modifies weak causality relation (eq. 1.11). Both phases 
result from M-matrix. 

In case of asymmetric barrier systems we can introduce times separately for reflection and 
transmission: 

, _ d((fi2+dk) . , _ d{ip2+dk+A(p) /-, nj\ 

^Monodromy, re f I — ) ^ M onodromy ,trns \ ') 



2 Nimtz experiments in view of the monodromy matrix 

The barriers in Nimtz experiment [12] consist of two photonic lattices which arc separated by 
an air gap. Each lattice consists between one and four equidistant Perspex layers separated by 
an air. The refractive index of Perspex is n=1.61 in the measured frequency region. In order to 
build a photonic barrier for the microwave signals, the thickness of the Perspex b=6.0(or 5.0) 
mm and the air layers a=12.0, (8.5) mm present a quarter of the microwave carrier's wavelength 
in barrier = c/(n/c) = 22.1, (20.4)mm and in air Aq = c//c = 35.5, (32.8)mm respectively. 
The air space dcav = 130, (189)mm between the two lattices forms a cavity and extends the total 
length of the barrier. The resonance frequencies of the cavity can be calculated on the base of 
monodromy matrix and are in case of two setups 1097. Mhz (or 764. MHz); according to Nimtz 
fres = c/{2d) is (1153) or (794MHz). 

The calculated transmission and wave function phases (according to monodromy for symmetric 
photonic lattices the reflection phase equals the transmission phase) are displayed in Fig(7,8). 

In Fig 8 we marked three areas with anomalous dispersion. 

Nimtz assumes that the frequency spectrum of the microwave signal lies completely in the 
nonresonant "forbidden" frequency region between 11 fres and 12 fres- Using the monodromy 
calculation method and if k^^^ is correct, fc should be shifted in comparison to Nimtz data (see 
figs 7,8) i.e. k from the value 0.1769 up to ~ 0.19 equivalent , fc' — 9GHz . The supcrluminal 
k-regions weakly depend on small k^^^ changes. However k^^^ should be determine from the 



internal Perspex structure. Sometimes for one barrier there is the sharp change in both ipi,(p2 
phases behaviour at A; = k^'^^ if barrier width is comparable or bigger then particle wave length . 
for 

Two barriers Nimtz experiment (9.15GHz) 
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FIGURE 7. Two barriers (9.15GHz) experiment: thickness of barrier is 5.0mm, dcav = 189mm. 
The superluminal speed changes gradually from 6.9c ai k = 0.12 to 3.4c near /c, cf. figS. Total 
width of the system is 199mm. 

Eight barriers (9.15GHz) 
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FIGURE 8. Eight barriers (9.15GHz) experiment at 9.15GHz: thickness of barrier is 5.0mm, 
of air layers is 8.5mm, while dcav = 189mm. The superluminal speed in region / is ~ 23c while 
in region II ~ 14c. Total width of the system is 280mm. 

It is not easy, from the phase curves to say where is placed. should be found by any 
independent method. 

At the end we present phase shift analysis in case of the Kiang model with 10 (^-barriers [10,11]. 
In allowed bands the "quantum hurdler" is retarded. But it is not so simple in case of forbidden 
bands. If the particle is reflected, negative slope can be related with penetration depth and for 
the transmitted wave with advance speed. 
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FIGURE 9. The Kiang model with ten 5 barriers is a good example of ("superluminal") time 
advance effect in the successive k- forbidden regions. As in [10], we used for that model fiA^ = 5, 
where is barrier penetrability and is the interbarrier distance. The slope of the phase line 
characterizes the superluminal speed i.e. for the deepest band we have ~ 19c for the next bands 
this speed decreases and is adequately: 10c, 5.5c, 4c and 2.5c in the latest band. In allowed 
regions the slopes are positive and give retardations. The group time delay is much bigger at 
resonances then between them - the phase curve oscillates strongly. 



2.1 Final remarks 

We assume that d is system width and t^^^ = d/{c) is the "classical" time needed to travel 
through barriers. For symmetric device both waves reflected and transmitted received the same 
phase shift <f2 (modulo tt /2) . The monodromy phase shift analysis of reflection and transmission 
amplitudes rejects reflection from the front and treats both processes as occurring after the time 
Ttun = + T^^^ where r^^^ results from (p2 ■ From the M- matrix point of view reflection and 
transmission delay if defined as Ttun = {d + d(p2/dk)/v = 2dpen/v is positive and allow us to 
interprete dpen as penetration depth in case of reflection. But we are unable to say from which 
position particle is transmitted. In that case interpretation is not easy. May be the tunneling 
wave function is strongly repelled from forbidden bands resulting in negative phase derivative. It 
seems that the tunneling particle needs much shorter time to travel through barrier than in free 
space. Numerical calculations (with M-matrix) show that always dpen > according to weak 
Wigner causality. If 2dpen is actual distance seeing by the tunneling particle then 2dpen/Ttun 
is typical speed (v) in the matter, but d/rtun gives the advance speed. We stress once more 
the transfer matrix enable us to find /c-dependence of phases. The phase shift 992 alone is not 
monotonic and corresponds to phase shifts from S matrix. Analysis of amplitudes from equivalent 
S matrix cause problems. We must define the arrival time for reflected particle and "departure" 
time for transmitted one. The probabilistic interpretation of R and T suggests that reflection 
occurs at barrier front [14] and then anomalous dispersion brakes weak Wigner causality, t^^^ 
time must be separately defined for the reflected and transmitted waves, in consequence there 
is no interference between these waves. When we calculate <f2ik) in allowed bands in simple 
(5-barrier systems we recognize that it is monotonic function of k with nontypical resonance 
structure of another origin. The quantum hurdler is retarded as it can be seen from (5argT in 
the Kiang model (fig. 9) but forbidden bands push out the particle. This structure emerges from 
interference effects between both ( reflected and transmitted or incoming and reflected) elastic 
channels in "continuum" including tunneling. Possible superluminal area are seen in Fig. (5- 9). 
i.e. anomalous dispersion < but it is questionable if ipi = const or is negligible - (^1 is 
oscillating function of wave number k . We analyze anomalous dispersion (without absorption 
or dissipation) only from two channel interaction point of view. In usual propagation of light in 
refractive media we reject influence of reflected wave. 

Brillouin has written in his book [15] "it is impossible to think of refractive medium without 



dispersion" (and energy loss) so the questions - what velocity coincides with elastic tunneling 
as well what the nature of anomalous dispersion is - remain open. 

We thanks to Nimtz for kind scientific cooperation. For stimulating discussion we are grateful 
to A.Horzela, V.S.Olkhovsky, E.Recami, S.Maydaniuk. 
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Figure descriptions 

FIGURE 1. Inward,outward, tunneling etc. waves in scattering. During scattering only one wave 
(here) outward or inward is modified. There is no cross terms between inward and outward fluxes . 
'^tun is not incorporated in nor in It is not clear if the reflected wave ^'re/K- \S\ = 1) 

is equal to vj/o"*"". In shadow region the complete wave function must vanish [16] (^/i = 0), there 
is no place for "^^tun in V'- 

FIGURE 2. The monodromy problem as defined in [6] for one barrier. Mutual relations be- 
tween ingoing and outgoing (from left or right side) particle wave functions are displayed. (In 
the picture bars mean complex coniugations.) The transition from initial state to final one 

( ^ . ( ^ . L.. ( m • 



is given by the unimodular M matrix , , easy de- 

V J^n \ J out V J 

ducted from transitions as drawn in picture. 

FIGURE 3. Set of barriers as used in tunneling and transmission or reflection. 

FIGURE 4. Transmission through one barrier three units [mm] wide, k^'^^ is barrier height in k 

units [1/mm]. —kd represents maximal negative phase slope according to weak Wigner causality. 

Slope of = 'P2 + 'Pi — kd is connected with the group velocity in transmission through the 

barrier. {d/vdk){(p2) = St is time "delay". The wave length A at AC(0) is smaller than the barrier 

width. 

FIGURE 5. Transmission through one Nimtz barrier 6mm wide, k^^^ is barrier height in k 
units (here mm~^). As in fig.4 (—kd) represents maximal negative phase slope according to 

weak Wigner causality. The wave length A at is bigger than the barrier width so phase 
characteristic is dominated by </?i. The slopes of (pi{k) and 6argT phase curves are positive and 
give retardation. However ip2 — kd = 2r] suggests small speed advance. 

FIGURE 6. Transmission through two equal barriers each one length unit [mm] wide. The 
cavity diameter is one millimeter wide too. k^'^^ is barrier height in k units [1/mm] the same as 
for one barrier tunneling (cf. fig.4). 

FIGURE 7. Two barriers (9.15GHz) experiment: thickness of barrier is 5.0mm, dcav = 189mm. 
The superluminal speed changes gradually from 6.9c at A; = 0.12 to 3.4c near fc, cf. figS. Total 

width of the system is 199mm. 

FIGURE 8. Eight barriers (9.15GHz) experiment at 9.15GHz: thickness of barrier is 5.0mm, 
of air layers is 8.5mm, while dcav = 189mm. The superluminal speed in region 7 is 23c while 
in region II ~ 14c. Total width of the system is 280mm. 

FIGURE 9. The Kiang model with ten 5 barriers is a good example of ("superluminal") time 
advance effect in the successive k- forbidden regions. As in [10], we used for that model QA^ = 5, 
where Q is barrier penetrability and is the interbarrier distance. The slope of the phase line 
characterizes the superluminal speed i.e. for the deepest band we have ~ 19c for the next bands 
this speed decreases and is adequately: 10c, 5.5c, 4c and 2.5c in the latest band. In allowed 
regions the slopes arc positive and give retardations. The group time delay is much bigger at 
resonances then between them - the phase curve oscillates strongly. 



